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of  ^  (analysed  to  second  order)  and  •  To  illustrate  the  compact  nature  of 
individual  expressions,  the  (untruncated)  perturbation  formulae  due  to  are 

given. 

This  is  the  text  of  a  paper  that  was  presented  (on  9  August  1989,  as  Paper 
AAS  89-451)  at  the  1989  AAS/AIAA  Astrodynamics  Specialist  Conference,  held  in 
Stowe ,  Vermont.  The  paper  is  printed  here,  as  pages  3-24,  in  the  format  required 
for  publication  of  the  Conference  Proceedings  in  Ado.  Astronaut.  Soi.  (Vol.71, 
pp  1229-1250,  1990).  In  an  edited  form,  it  has  also  now  been  published  in  the 
Journal  of  the  Astronautical  Sciences  (Vol.39,  pp  65-85,  1991).  The  paper  was 
originally  prepared  as  a  greatly  shortened  version  of  RAE  Technical  Report  89022. 

Copyright 

© 

Controller  HMSO  London 
1991 


UNLIMITED 


LIST  OP  CONTENTS 


Page 

INTRODUCTION  3 

POTENTIAL  DEVELOPED  VIA  THE  INCLINATION  FUNCTIONS  5 

ECCENTRICITY  FUNCTIONS  REQUIRED  IN  SUBSEQUENT  ANALYSIS  7 

RATES  OF  CHANGE  OF  OSCULATING  AND  MEAN  ELEMENTS  9 

COORDINATE  PERTURBATIONS  (GENERAL  CASE)  1] 

The  Perturbation  6r  12 

The  perturbation  6b  13 

The  perturbation  6w  -14 

THE  SPECIAL  CASES,  AND  INTEGRATION  CONSTANTS  17 

Mandatory  Constants  for  6a  17 

Constants  for  6e  and  5M  17 

Constants  for  6i  and  6£2  18 

Forced  Terms  in  6w  18 

Constants  for  6cj  19 

RESULTS  FOR  J.  19 

4 

EXTENSION  TO  TESSERAL  HARMONICS  22 

CONCLUSION  23 

References  23 

Report  documentation  page  inside  back  cover 


3 


AAS  89-451 


PERTURBATIONS,  UNTRUNCATED  IN  ECCENTRICITY, 
FOR  AN  ORBIT  IN  AN  AXI-SYMMETRIC 
GRAVITATIONAL  FIELD 

R.  H.  Gooding 

By  recourse  to  a  particular  definition  of  a 
satellite’s  mean  orbital  elements  and  to  a 
particular  system  of  spherical  polar  coordinates 
based  on  the  mean  orbital  plane,  an  orbital  theory 
has  been  developed  that  leads  to  extremely  compact 
first-order  perturbation  formulae  associated  with 
the  general  zonal  harmonic,  .  The  formulae  are 
complete  (untruncated  in  eccentricity)  and 
generalize,  via  recurrence  relations,  the  author's 
earlier  results  for  the  effects  of  J2  (analysed  to 
second  order)  and  J3  .  To  illustrate  the  compact 
nature  of  individual  expressions,  the  (untruncated), 
perturbation  formulae  due  to  Jjj  are  given. 

INTRODUCTION 

In  an  earlier  paper1,  the  author  summarized  a  theory  for  satellite 
perturbations  due  to  the  harmonics  J2  and  J3  of  the  Earth's  gravit¬ 
ational  field,  presenting  formulae  that  are  complete  (untruncated  in 
eccentricity)  to  second  order  in  J2  and  first  order  in  J3  .  The 
novelty  of  the  theory  arises  from  the  way  in  which  short-period  pertur¬ 
bations  in  the  osculating  element  are  amalgamated  into  perturbations  in 
a  set  of  spherical- polar  coordinates  (r,  b,  w),  based  on  a  mean  orbital 
plane:  r  is  the  geocentric  distance  (radial  direction),  whilst  b 
and  w  are  quasi-latitude  (cross-track  direction)  and  quasi-longitude 
(along-track  direction).  Particular  definitions  of  the  mean  orbital 
elements  were  adopted,  to  make  the  coordinate-perturbation  expressions 
as  compact  as  possible,  these  expressions  being  complemented  by 
formulae  for  the  rates  of  change  of  the  mean  elements,  to  take  care  of 
the  secular  and  long-period  behaviour. 

Full  details  of  the  J2/J3  theory  are  given  in  a  recent  RAE 
report2,  and  a  subsequent  report^  gives  the  details  of  the  theory’s 
extension  (to  first  order)  to  the  general  zonal  harmonic,  .  The 
present  paper  is  essentially  a  precis  of  Ref.  3,  which  will  often  be 
referred  to  simply  as  'the  Report'. 
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The  starting  point  for  the  theory  is  the  decomposition  of  Uj  , 
the  potential  due  to  Jj,  »  into  a  finite  sum  of  component  terms,  uj£  , 
each  of  which  is  associated  with  a  particular  function  of  the  orbital 
inclination.  Then,  to  integrate  Lagrange’s  planetary  equations, 

(p/r)^-1  (where  p  is  the  semi-latus  rectum  of  the  orbit)  is  expanded 
as  a  finite  sum  of  terms,  each  of  which  is  associated  with  an 
eccentricity  function.  The  two  sets  of  functions  can  be  generated  by 
the  use  of  recurrence  relations,  and  certain  of  these  relations  are 
also  vital  in  the  development  of  the  coordinate-perturbation  formulae 
(expressions  for  6r,  6b  and  6w,  associated  with  u£  )  that  constitute 
the  principal  results  of  the  paper. 

For  each  of  6r,  6b  and  6w,  the  general  formula  (associated  with 
)  involves  a  sum  of  terms  over  a  third  index,  j  .  Certain  values 
of  j  lead  to  zero  denominators  in  particular  terms,  but  for  6r  and 
6b  these  terms  can  be  eliminated  altogether,  by  the  definition  of  the 
mean  elements.  For  6w  ,  on  the  other  hand,  most  of  the  potentially 
infinite  terms  do  not  just  disappear;  instead,  they  are  replaced  by 
particular  finite  terms  that  are  induced  by  the  element  definitions. 

The  formulae  for  these  'particular'  terms  supplement  the  general  formu¬ 
lae.  A  further  complication  originates  from  the  necessity  to  change 
the  integration  variable,  in  the  planetary  equations,  from  t  (time)  to 
v  (true  anomaly).  Since  the  mean  elements  evolve  with  t  ,  rather  than 
v  ,  an  additional  short-period  perturbation  is  induced  by  every  secular 
or  long-period  component  of  the  rate  of  change  of  a  mean  element;  the 
perturbations  induced  by  the  long-period  components  (but  not  the 
secular  components)  are  best  treated  like  the  pure  short-period  pert¬ 
urbations,  by  amalgamating  them  into  perturbations  in  the  coordinates. 

To  exemplify  the  general  theory,  specific  results  are  given  for 
the  harmonic  .  The  secular  and  long-period  expressions  for  the 
element  rates  of  change  are  well-known,  but  the  expressions  for  6r  , 

6b  and  6w  have  not  been  given  before. 

The  formulae  associated  with  the  general  zonal  harmonic  do  not 
immediately  extend  to  the  tesseral  harmonics,  because  of  the  complica¬ 
tion  introduced  by  the  Earth's  rotation.  For  a  zero  rate  of  rotation, 
however,  the  extension  is  simple  and  straightforward,  if  of  little 
practical  value,  and  this  is  the  final  topic  of  the  paper. 

It  is  convenient  to  conclude  the  introduction  with  some  remarks  on 
notation.  We  assume  the  usual  elliptic  elements  (a,  e,  i,  fi,  u,  M), 
with  n  and  o  defined  such  that  M  -  o  +  Jn  dt  ;  further,  this  inte¬ 
gral  will  often  be  written  just  as  /  .  As  in  the  author’s  earliest 
work*4  on  orbital  theory,  we  also  find  it  convenient  to  utilize  the 
quasi-elements  ,  p  and  L  ,  defined  such  that  dxp  -  du>  +  c  dfl  , 
where  c  -  cos  i  (and  we  also  write  s  =  sin  i  ),  dp  -  do  ♦  q  dij» 
(where  q2  «  1  -  e2)  and  dL  -  dM  *  q  dip  .  Finally,  and  as  a  valuable 
shorthand,  we  define 

Cj  -  cos  (jv  +  ku')  and  Sj  -  sin  (jv  +  ku’)  ,  (1) 
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where  u'  -  oa  +  v  -  -yif  (argument  of  latitude  measured  from  the  north 
apex  of  the  orbit);  when  no  ambiguity  is  possible,  we  will  frequently 
omit  the  superfix  k  in  this  notation. 

POTENTIAL  DEVELOPED  VIA  THE  INCLINATION  FUNCTIONS 


The  standard  expression  for  U£  is 

U£  =  -  ^  J£  (R/r)*  P£(sin  6)  . 


(2) 


As  in  Ref.  we  expand  P^sin  B)  via  the  addition  theorem  for  zonal 
harmonics;  thus 


P£(sin  g)  =  l  uk  ((*  ~  k-jj  Pj[(0)  Pj(c)  cos  ku* 


(3) 


Here  uq  =  1  •  uk  =  2  if  k  >  0  ,  and  the  Legendre  function  P£  is 
defined  by 


P  £  ( c )  -  s^ 


dkP£(c) 

dck 


(*0 


The  second  factor  (the  k'th  derivative)  in  (4)  is  a  polynomial  in 
c  ,  which  (with  k  i  ()  does  not  vanish  when  c  ■>  1  ,  its  value  then 
being  ( i.  +  k)!/{2k  k!  (£,  -  k ) ! }  .  Hence  this  factor  may  be  normalized, 
in  a  certain  useful  sense,  and  we  write 


dkP£(c)  (£,  +  k)!  .k/ ,  \ 

dck  2k  k!  (£  -  k)!  l 


(5) 


where  A^(i)  is  a  pure  polynomial  in  s2  if  k  has  the  same  parity 
as  l  ,  but  has  an  additional  factor  c  if  k  and  i,  are  of  opposite 
parity;  in  each  case  the  constant  term  in  the  polynomial  is  unity,  by 
the  normalization.  Explicit  expressions  for  the  A£(i)  are  given  in 
Table  1  of  the  Report. 

We  can  now  rewrite  (3)  as 


l 

P£(sin  6)  =  l  alk  sk  A^( i )  cos  ku’  ,  (6) 

k-0 

where  the  constant,  a£k  ,  is  given  by 

°£k  “  uk  /<2k  k!)  .  (7) 

It  is  clear,  from  the  last  paragraph,  that  P£(0)  vanishes  when  k 
and  i  are  of  opposite  parity,  and  it  may  be  shown  that  when  the 
parity  is  the  same, 
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P$(0) 


t(  l-k ) 


_ U  *  k)! _ 

2*  {*(£  ♦  k ) }  I  {*(£  -  k)}! 


In  substituting  (6)  into  (2)  it  is  of  great  benefit  to  introduce  a  new 
quantity,  A^k  ,  defined  by 


(R/p)£  ci«|,  sk  Ap  (i )  . 


Tnis  permits  us  to  write 


u£  -  -  |  (p/r)*+1  A£k  C&  .  (1( 

1/  .  . 

It  will  be  noted  that,  whereas  A^(i )  is  defined  and  useful 
regardless  of  parity,  Aik  (and  hence  )  is  only  non-zero  when  k 
and  2.  are  of  the  same  parity.  However,  a  use  will  be  found  for  quan¬ 
tities  that  behave  in  the  opposite  way  from  a £k  and  A^k  ,  antici¬ 
pating  which  we  define  (with  bold  letters  to  make  the  distinction) 

al<  =  U  -  k  ♦  1)  p£+1(0)/(2k  <1  £)  (11 


A£k  =  J£  (R/p)*  alK  s<  a£(  i )  , 


where  k  has  been  replaced  by  k  to  signify  that  we  now  have 
quantities  that  are  non-zero  only  when  k  and  £  are  of  opposite 
parity. 

We  will  require  derivatives  of  the  inclination  functions.  It  is 
evident  from  (5)  that 


aruiu» 


U  -  k)(£  +  k  +  1)  _  .k-n,,s 
- 2(k  +  1 ) -  5  A*  (1)  ; 


from  this  and  (9)  it  follows  that  the  (partial)  derivative  of  A^k 
with  respect  to  i  is  given  by 

A£k  *  J£  (R/P)*  a£k  {kc  A£(i)  -  U  ~  k2(kViV  1 }  s2  Ar1(i)}  •  (1J° 

We  will  also  require,  finally,  the  particular  combinations  of  A^k  and 
A^k  denoted  by  A^k  and  Ajk  ,  and  given  by 

A^k  m  ks  ^  A^k  ±  c  ^  A^k  i  (1 5) 

these  definitions  reduce^  to 

4jk  -  (t  ~  ^’Vi)  ***’  «$*'»>  <’6> 
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Ajk  =  2k  J£(R/p)*  a£k  c_1  sk_1  A^_ 1 C i )  .  (17) 

Various  recurrence  relations  for  the  functions  A^(i)  were  given 
in  the  Report,  and  these  may  be  used  (if  any  possibility  of  instabil- 
ity5  is  disregarded)  to  generate  the  A£k  and  k^K  .  It  is  also 
possible  to  connect  A^k  and  A£k  with  the  appropriate  k^K  ,  and  the 
following  pair  of  formulae  will  be  needed  in  the  sequel: 


and 


Aft,  k  +  1  _  Aft,k-1 

uk+1  uk-1  • 


2kcs  1 


Aft,  k  +  1 
uk  +  1 


+  Aft, k- 1 

uk-1 


(18) 


(19) 


ECCENTRICITY  FUNCTIONS  REQUIRED  IN  SUBSEQUENT  ANALYSIS 

The  term  Uj,  of  the  potential,  specified  by  (2),  has  now  been 
decomposed  into  the  defined  by  (10),  the  latitude  ( B)  having  been 

eliminated.  The  longitude  was  absent  from  U£  from  the  beginning, 
because  of  axial  symmetry,  so  it  remains  to  eliminate  the  radius  vector 
(r)  .  Since  p/r  =  1  +  e  cos  v  ,  an  expansion  of  the  form 

ft-T 

(p/r)^-1  =  I  uj  B£j  cos  Jv  (20) 

j=0 

is  possible,  for  ft  2  1  ,  and  we  regard  B£,  as  defined  by  this  expan¬ 
sion;  clearly,  B£,  is  a  polynomial  in  e  .  We  shall  find  it  useful, 
and  entirely  natural,  to  extend  the  definition  of  B£j  to  negative 
j  ,  by  defining  Bjj  »  Bft,|j|  »  and  to  take  B$j  =  0  when  I J I  *  A  • 

On  this  basis  we  can  replace  (20)  by 

(p/r)*'1  -  l  B£j  C<?  ,  (21) 

where  the  summation  effectively  runs  from  j  =  -®  to  j  -  +•  .  The 
B £ ^  are  directly  related  to  the  Hansen  X  functions  of  classical 
celestial  mechanics,  since 

Bjj  =  q2*"1  Xq*'1^  .  (22) 

We  proceed,  as  in  Ref.  H,  to  express  the  e-polynomial  B^  ,  when 
ft  >  1  and  0  <  j  <  ft  ,  in  terms  of  a  normalized  polynomial,  tne 
connecting  relation  being 

Bftj  *  (* j  ’)  (e/2) J  B{(e)  ;  (23) 

B|(e)  in  (23)  is  a  polynomial  in  e2  ,  with  constant  term  unity  by  the 
normalization.  Explicit  expressions  for  the  B|(e)  are  given  in 
Table  3  of  the  Report. 
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In  contrast  with  the  A£k  ,  however,  it  is  usually  better  to  work 
with  the  By  directly,  rather  than  through  B|(e)  and  (23).  One 
reason  for  this  is  that  only  alternate  values  of  the  A£k  are  non¬ 
zero,  whereas  (for  |  j  |  <  £  and  e  *  0)  all  the  By  are  non-zero. 
Further,  no  difficulty  arises  with  the  By  when  j  <  0  ,  whereas 
B|(e)  would  then  be  infinite  (if  J  j |  <  £).  We  can  even  allow  £  to 
be  negative  (or  zero)  as  well  as  j  ;  the  validity  of  this  follows 
from  the  universality  of  (22). 

In  regard  to  derivatives  of  the  eccentricity  functions,  it  can  be 
shown  that,  for  1  $  j  <  2.  , 

^-jBj(e)}  =  2je_1  {BJy](e)  -  B*[(e)}  .  (2H) 

The  universal  formula  (valid  for  all  j  )  for  the  derivative  of  By 
can  then  be  shown  to  be 

Bij  -  H-IJ-I  -  J  e''  B£j  •  (25> 

However,  because  we  only  introduce  the  By  after  each  planetary 
equation  has  been  set  up,  we  effectively  only  use  (25)  in  expressing 
the  rates  of  change  of  the  mean  elements.  Since  this  involves 

■h  <q  4IK  Blk>  ■  q'1  Hk  <12  Bik  *  (2*  -  ')  «  Blk)  •  <«> 

we  define 

Elk  “  q2  B£k  ♦  (21  -  1)  e  B£k  ;  (27) 

then  (25)  leads  (via  a  recurrence  relation)  to 

E£k  =  e-1  (£e2  -  k)  B£k  +  (£  -  k)  B4fk_1  .  (28) 

A  number  of  recurrence  relations  for  the  By  are  given  in  the 
Report.  However,  these  may  all  be  generated  fromJjust  two  relations, 
and  we  use  here  the  following  pair,  as  being  most  useful  in  the  sequel: 

(j  -  £)  e  B£y_1  ♦  2j  By  ♦  (j  +  £)  e  Byj  +  1  -  0  (29) 

and 

( £  ~  j)  ©  B£+i  , j  +  £q2  ®£,J  +  1  ~  +  J  +  ^)  ®£+1  , j+1  *  0  .  (30) 

Alternatively,  we  can  replace  the  unsymmetrical  (30)  by  the  more 
symmetrical  (and  simpler)  equation 

ie  (B£,j-1  "  B£,j  +  1)  “  2J  B£+1,j  •  (31) 

given  by  Zaf iropoulos^  in  his  orbital  theory  for  an  axi-symmetric 
field. 


TM  Sp  379 


9 


RATES  OF  CHANGE  OF  OSCULATING  AND  MEAN  ELEMENTS 

It  is  straightforward  to  derive  finite  expressions  for  the 
perturbations  in  the  osculating  elements,  so  long  as  we  change  the 
variable  in  Lagrange's  planetary  equations  from  t  to  v  before 
integrating.  The  equation  needed  for  the  change  of  variable  is 

dv/dt  -  n  q“3  (p/r)2  .  (32) 

We  obtain  the  expression  for  6a  first,  however,  because  a  direct  deri¬ 
vation  is  available,  which  is  simpler  than  using  the  planetary  equation. 

The  starting  point  for  the  direct  derivation  of  6a  is,  as  in 
previous  work1-1*,  the  existence  of  an  energy-related  absolute  constant 
of  the  motion,  which  we  denote  by  a'  ;  the  relation 

a  =  a'  (1  +  2aU/u)  (33) 

is  exact  for  any  time-independent  disturbing  function,  U  ,  and  in 
particular  for  the  axi-symmetric  .  It  follows  that  there  is  no 
long-term  variation  in  a  ,  to  whatever  order  of  magnitude  the  pertur¬ 
bation  analysis  is  conducted.  Further,  the  short-period  perturbation, 

6a  ,  is  given  exactly,  on  substituting  for  u£  from  (10);  thus 

6a  -  -  2a'q-2  A^k  (p /r)l*]  Cq  .  (3*0 

This  does  not  mean  that  an  exact  perturbation  can  be  written  down  for 
semi-major  axis,  however,  as  the  right-hand  side  of  (3^)  is  expressed  in 
terms  of  osculating  elements;  as  soon  as  mean  elements  are  introduced, 
the  result  is  no  more  than  a  first-order  perturbation  expression. 

To  present  6a  in  the  form  appropriate  for  later  use,  we  combine 
Cq  with  one  of  the  factors  p/r  .  We  retain  another  p/r  factor  expli¬ 
citly,  and  expand  the  remaining  (p/r)^-1  by  (21).  By  this  means,  we 
effectively  transform  the  term  2Cq  ,  for  example,  into  Cj  +  C_j  .  But 
each  pair  of  terms  (such  as  this)  for  positive  j  ,  in  the  infinite 
summation  of  (21),  is  matched  by  the  same  pair  (in  reverse  order)  for 
negative  j  ,  so  we  can  express  our  desired  result  as 

6a  **  -  aq-2  A£k  (p/r)  l  Bjj  (eC i  +  2Cj  +  eCj+1)  .  (35) 

For  the  remaining  elements  (e,  i,  ft,  w  and  M)  ,  it  is  convenient 
to  start  with  the  time  rates  of  change  of  p,  pc2,  fl,  iji  and  p,  since  the 
right-hand  sides  of  the  planetary  equations  for  these  rates  are  the 
single-term  expressions  (2q/na)  3U/3aj,  (2qc/na)  3U/3fl,  (1/na2qs)  3U/3i, 
(q/na2e)  3U/3e  and  (-2/na)  3U/3a,  respectively.  We  now  use  (32)  to 
change  the  variable  to  v  ,  and  (21)  to  eliminate  r  ,  as  a  result  of 


which  we  get  (with  only 
infinite  summation): 

a  finite  number  of  non-zero  terms  from  each 

dp/dv 

-  2kp  Aik 

1  Sq  » 

(36) 

d(pc2)/dv 

-  0  , 

(37) 
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dfi/dv 


(38) 


di^/dv 


s  1  A£k  E  B£j  Cj  > 

~  V  e_1  A£k  l  By  1  ( «-  +  1  -  k)e  Cj_2  ♦  2(1  *  1  -  2k)Cj_1 
+  2 (£  +  1  )e  Cj  +  2(£  +  1  ♦  2k)Cj  +  1  ♦  (£  +  1  +  k)e  Cj+2}  (39) 


and 

dp/dv  =  -  2(1  +  1)q  A£k  l  By  Cj  .  ( 40) 

We  also  require  a  formula  for  the  rate  of  change  of  the  (shorthand) 
quantity  /  ;  the  appropriate  formula  is 


d//dv  =  n’/V  ♦  3q  A£k  l  By  Cj  , 
where  n'  is  the  exact  constant  defined  by  n'2  a’3  =  y  . 


(^1) 


To  obtain  the  rates  of  change  of  the  mean  elements,  we  select  the 
terms  involving  C^k  or  S^k  (in  the  equations  other  than  (39)  this 
implies  just  j  -  -k  ) ,  since  these  are  the  terms  that  are  independent 
of  v  .  We  abandon  p  and  pc2  in  favour  of  e  and  i  ,  but  we 
retain  ij*  and  p  ,  since  the  secular  and  long-period  rates  of  these 
(mean)  ’elements'  are  useful  in  practice,  and  this  is  also  true  for  L  ; 
since  ~5  (i.e.  a'  )  does  not  vary,  our  full  results  may  be  written 
(after  simplification  where  necessary): 


^£k 

=  -  kne- 

1q2  A£k  B£k  S-k  * 

(42) 

T£k 

■  kncs-1 

A£k  B£k  s-k  * 

(43) 

^£k 

=  -  ns"1 

A£k  B£k  c-k  » 

(44) 

*£k 

=  -  ne-1 

A£k  E£k  c-k  » 

(45) 

^£k 

-  -  2(£ 

+  1)nq  A£k  B£k  C-k 

(46) 

L£k 

=  -  (2£ 

-  D  nq  A£(<  B£k  C*k  . 

(47) 

The  remaining  terms  from  (36)  -  (*10),  those  that  depend  on  v  , 
lead  to  the  pure  short-period  pertur bations  in  the  elements.  For  future 
reference  we  supplement  this  set  of  equations  by  giving  those  for  e  , 
i  ,  M  and  L  ,  which  can  be  derived  from  the  original  set  (with  an 
appropriate  formula  for  da/dv  included);  thus: 


de/dv 


-J-  A£k  l  By  { ( k  -  £  -  1  )e  Sj_2  +  2(2k  -  £  -  1 ) S j _ , 

+  6ke  Sj  ♦  2(2k  +  £  +  1)Sj  +  1  +  (k  ♦  £  +  i )e  Sj+2)  . 

kcs-1  A£k  l  By  Sj  , 


(48) 


di/dv 


(49) 


dM/dv 


(50) 


1  1 

=  -q-  e_1q  A£k  l  B^j  {(£  +  1  -  k)e  Cj_2  +  2(1  ♦  1  -  2k)Cj_1 
-  6(2,  -  1)e  Cj  +  2(2.  +  1  +  2k)Cj+1  +  (2  +  1  +  k)e  Cj  +  2} 

and 

dL/dv  =  -  (22  -  1)q  A^k  l  Bjj  Cj  .  (51) 

The  integration  of  these  equations  is  entirely  straightforward,  since 
the  analysis  is  only  being  taken  to  first  order;  also,  zero  denominators 
cannot  occur,  because  the  terms  that  would  produce  them  are  precisely 
the  ones  that  have  been  dealt  with  separately  and  are  no  longer  present. 

To  save  space,  therefore,  these  integrals  are  omitted. 

In  the  rest  of  the  paper  we  are  concerned  with  the  amalgamation  of 
the  integrals  for  short-period  6a  etc  into  formulae  for  the  coordinate 
perturbations ,  6r,  6b  and  6w.  However,  we  conclude  the  present  section 
with  two  remarks  about  the  secular  and  long-period  perturbations  -  they 
will  be  relevant  in  obtaining  the  formulae  associated  with  ,  to  ex¬ 
emplify  the  general  results.  First,  the  long-period  perturbations 
induce  additional  terms  in  6r,  6b  and  6w,  as  noted  in  the  Introduction. 
Second,  when  2  is  even  it  is  convenient  to  deal  with  the  secular  pei — 
turbation  in  M  by  use  of  a  'mean  mean  motion' ,  n,  that  is  not  the 
same  as  n'  (except,  as  it  happens,  for  2  =  2);  since  ~S  is  ident¬ 
ified  with  a'  ,  this  means  that  Kepler's  third  law  does  not  hold  for 
mean  n  and  a  .  (This  is  covered  in  detail  in  the  Report,  and  Ref.  7 
is  almost  entirely  devoted  to  the  matter.) 

COORDINATE  PERTURBATIONS  (GENERAL  CASE) 

In  this  section  we  develop  general  expressions  for  the  6r  ,  6b 
and  6w  that  can  be  derived  from  the  (untruncated)  first-order  formulae 

6r  =  (r/a)  6a  -  (a  cos  v)  6e  +  (aeq-1  sin  v)  6M  ,  (52) 

6b  =  (cos  u')  6i  +  (s  sin  u')  6P  (53) 

and 

6w  =  5i|>  +  {q-2  sin  v  (1  +  p/r)}  6e  +  q~3  (p/r)2  6M  .  (5*0 

Special  cases,  associated  with  the  particular  choices  of  mean  elements 
that  will  be  derived,  are  reserved  for  the  next  section  of  the  paper. 

Generation  of  the  expressions  for  6r  and  6w  is  essentially 
straightf orward  in  that  the  analysis  starts  with  the  6a  etc  due  to 
Uy  and  finishes  with  6r£k  and  6w^k  .  With  6b  ,  however,  there  is  a 
complication,  due  to  the  appearance  of  u'  in  (53),  as  opposed  to  v 
in  (52)  and  (5*0;  we  deal  with  the  difficulty  by  deriving  6by_K  ,  rather 
than  6b£k  ,  where  <  has  values  of  opposite  parity  to  those  of  k  . 

We  do  not  give  expressions  for  6r,  6b  and  6w,  but  they  are  immedi¬ 
ately  available2* 3  from  the  expressions  for  6r,  6b  and  6w,  just  by 
replacing  Sj  and  Cj  by  (k  +  j)  n  Cj  and  -(k  +  j)  n  Sj  . 
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We  can  do  better  than  this  if  we  allow  for  the  (overall)  rate  of  change 
of  u5  ,  replacing  (k  +  j)  n  by  (k  +  j)  n  +  k  ,  assuming  Cj  and 
Sj  still  to  be  shorthand  for  Cj  and  S*j  . 

The  Perturbation  6r 

We  have  to  apply  (52)  with  6a  ,  5e  and  6M  given  by  (35)  and 
(the  integrals  of)  (48)  and  (50).  The  integrals  combine  in  a  very 
natural  way,  and  we  are  able  to  write 


6r 


"  *  a  ASlk  l  B 


Rj 


(55) 


where 

Rh 


=  e(j  l  -  1) 


[k  +  j  -  2 

.  *1  -  1  '  1 1c,  * 


k  +  j  +  1 


k  +  j 
e(j  -  l 


■'j-l 


+  2 


2j  +  &  -  1 
k  +  j  -  1 


+  U[k  J  j  +  k  +  j  +  2  jCj  +  1  * 


(56) 


It  can  be  seen  that  (55)  is  a  summation  in  which  Rj  ,  as  given  by 
(56),  has  three  components,  each  component  being  expressed  as  the  sum  of 
two  multiples  of  the  same  'C  quantity'.  We  separate  the  first  multiple 
from  the  second  (in  each  component  of  R j  ) ,  feeding  them  back  separate¬ 
ly  into  the  summation  of  (55);  this  leads  to  tv  istinct  summations 
that  we  can  denote  by  E_  and  1+  .  Thus  I_  involves  Z  By  R^_  ,  where 


=: 


j  ♦  l  -  1 
k  +  j  -  2 


eC 


j-1 


+ 


2 


2J  +  i.  -  1  r 
k  +  j  -  1 


J  -  a  +  i 


eC-i  + 


j  +  1 


(57) 


Since  sums  over  By  can  be  regarded  as  running  from  -»  to  +®  , 
it  follows  that  we  can  rearrange  the  three  sets  of  terms  in  ZBy  Rj_ 
so  that  (with  j  now  used  in  a  different  way) 

l  By  Rj_  -  H(j  +  a)e  BlJ  +  1  +  2(2j  ♦  £  -  1)By 

♦  3( j  -  £)e  Btjj_1}(k  +  j  -  D-1  Cj  .  (58) 

Using  the  recurrence  relations  (29)  and  (30),  we  can  simplify  this  to 

l  By  Rj_  -  2U  "  1)q2  I  (k  +  j  -  D'1  By1(j  Cj  .  (59) 

Similarly, 

I  By  Rj+  =  -  2(1  -  1)q2  l  (k  ♦  j  ♦  I)”1  Byyj  Cj  .  (60) 

The  final  result  we  require  now  follows  from  (55),  (59)  and  (60). 
Because  of  its  importance,  we  write  Cj  in  full.  Thus 

«rtk  •  -U  -  D  P  Alk  i  (k  .  j  .  .  j  -  ,)  =°s  <k“’  *  J»>  • 

j  ....  (61) 
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If  j  =  -k  ±  1  ,  there  is  a  zero  denominator  in  (61),  and  terms 
with  these  values  of  j  must  be  excluded;  in  the  next  section  we 
determine  integration  constants  for  6eik  and  6M£k  so  that  the  terms 
with  these  two  values  of  j  are  forced  to  zero.  It  will  be  noted  that 
all  the  cosine  terms  occurring  in  (61),  for  a  given  and  all  possible 

k  ,  are  distinct,  except  that  if  k  =  0  (£  even)  then  equal  and  opposite 

values  of  j  lead  to  identical  terms  in  cos  jv  .  On  this  basis  we  can 
derive  the  total  number  of  terms  required  to  express  6r  for  a  given 
value  of  £  .  This  is  done  in  detail  in  the  Report,  the  result  being 
£2  -  -£-(3£  -  1 )  when  £  is  odd  and  £^  -  ■?( 3£  -  2)  when  £  is  even. 

The  perturbation  6b 

We  get  6b  from  (53).  where  6i  and  6Q  are  given  by  the 

integrals  of  (A9)  and  (38).  This  is  on  the  assumption  that  6b  (=  6bjik) 

is  associated  with  ,  following  the  decomposition  of  Uj,  .  We  shall 
find,  however,  that  it  is  much  more  convenient  to  decompose  the  total 
6b  (associated  with  )  as  l  6b ,  where  the  summation  is  for  values 
of  tc  that  are  of  opposite  parity  to  £  and  we  no  longer  associate  the 
individual  5b  (=  6b£<)  with  specific  components  of  . 

In  relation  to  ,  we  get 

6b£k  -  -  l  B£j  {ir?-J  A£k  Cj  cos  u’  +  FT“J  A£k  sj  sin  u'}  *  (62> 

The  trigonometrical  products  are  replaced  by  sums,  in  the  usual  way,  and 
we  can  then  invoke  the  notation  of  (15)  to  write 

6b£k  -  -  4c  l  (k  ♦  jVUA$k  Of  ♦  A*k  C$+1>  .  (63) 

This  expression  may  be  contrasted  with  (55)  and  (56)  for  6r  .  In  view 
of  the  difference  in  superfix  (whereas  the  suffix  alone  varied  in  the 
terms  of  R j  ) ,  in  the  two  C  terms  of  (63),  we  would  now  like  to 
combine  a  pair  of  terms  with  different  k  indices,  before  the  summation 
over  the  j  index  operates. 

With  the  philosophy  just  indicated,  we  make  the  new  decomposition 

6b£  a  I  6b£*  -  (6A) 

where  each  6b  (=  6b1|C  )  is  of  the  form 

6b  -  l  Tj  Bgj  Cj  (65) 

and  we  require  an  expression  for  Tj  .  We  note  first  that  since  (for 
non-trivial  results)  k  runs  from  0  or  1  to  £  (taking  alternate 
values),  it  follows  that,  in  principle,  k  runs  from  -1  or  0  to 
£  +  1  (again  alternate  values,  but  of  opposite  parity  to  k):  for  the 
minimum  value  of  k  ,  only  the  term  in  A^k  ,  in  (63),  contributes  to 
Tj  ,  whilst  for  the  maximum  value  of  k  ,  only  the  term  in  A£k  contri¬ 
butes;  for  intermediate  values  (if  any),  both  terms  contribute.  But  we 
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can  straight  away  dismiss  the  'maximum  value'  (<  =  £  +  1),  because  Aj^ 
is  just  a  multiple  of  s*1  ;  from  this  it  follows  that  A^  ,  defined  by 
(15),  is  zero.  (Also  P^  =  0  anyway!)  We  find  that  we  do  not  require 
the  'minimum  value'  (<  -  -1)  either.  Then,  using  (16)  and  (17)  for  A£k 
and  A£k  ,  respectively,  we  get 

T  f  <  +  1  ( £  -  k  +  1 ) (£  +  k)  „ 

Tj  ■  J£  +  j  ~  1  ai,<+1  +  4tc(<  +  j  -  1)  £><-1 

x  (R/p)5,  sK  A^(i)  .  (66) 


The  quantity  in  curly  brackets  in  (66)  is  a  pure  constant,  in 
which  the  alk  are  given  by  (7):  thus  the  first  a  involves  p£+1 (0) 
and  the  second  involves  P£'  (0)  ,  these  being  given  by  (8).  By  rela¬ 
ting  these  to  p£+1(0)  ,  and  hence,  via  a£|C  given  by  (11),  to  A£(C 
given  by  (12),  we  can  rewrite  (66)  as 


*  *£> 
2u*. 


u<+i 


u<-1 


<  +  J 


K  +  j  "  1 


(67) 


The  preceding  argument  is  the  general  one,  for  1  S  k  £  £  -  1  , 
and  for  <  >  2  we  can  obviously  cancel  uK  out  with  uK+i  and  u^-i  . 
Very  conveniently  (with  the  full  argument  given  in  the  Report),  (67)  is 
also  correct  without  the  occurrences  of  u  for  <  »  0  ( £  odd),  and  even 
for  <  =  1  ( £  even)  on  the  basis  that  the  case  <  «  -1  is  then  auto¬ 
matically  covered. 


We  can  now  write  down  the  final  result  we  require,  on  substituting 
(67)  into  (65)  and  expressing  in  full.  Thus 

»u  ■  - »  Hk  i  u  *  j  .  n  Btj 003  UlJ'  * Jv)  •  <68) 

j 

As  already  indicated,  this  formula  is  unlike  (61),  the  correspond¬ 
ing  one  for  6r  ,  in  that  it  cannot  be  taken  in  isolation  as  relating  to 
a  sub-component  of  .  It  is  like  (61)  in  one  respect,  however,  in 
that  terms  of  Sbj^  with  j  »  -k  ±  1  are  excluded.  In  the  next  sec¬ 
tion  we  determine  constants  for  6iik  and  6«£k  (k  ,  not  k  ,  now 
being  the  appropriate  symbol)  such  that  these  terms  are  forced  to  zero. 

On  this  basis  we  can  derive  the  total  number  of  terms  (without  duplica¬ 
tion  of  Cj  )  required  to  express  6b  for  a  given  value  of  £  .  The 
result  (derived  in  the  Report)  is  £2  -  •!-(£  -  1)  when  £  is  odd  and 
£2  -  i(3£  -  2)  when  £  is  even  (the  latter  being  the  same  as  for  6r  ). 

The  perturbation  6w 

The  analysis  for  6w  is  much  more  like  the  6r  analysis  than  the 
6b  analysis,  because  each  Uq  can  again  be  treated  separately  through¬ 
out.  There  are  two  complications,  however.  First,  (54)  effectively 
involves  cos  2v  and  sin  2v  ,  not  Just  cos  v  and  sin  v  (we  see  this 
at  equation  (69),  following),  and  this  means  that  the  values  j  -  -k  ±  2 
are  special  as  well  as  j  -  -k  ±  1  .  Second,  we  cannot  take  6w  to  be 
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zero  for  any  of  these  special  cases,  since  the  constants  in  6e  and  6M 
must  now  be  assumed  to  have  been  already  assigned;  formulae  for  the  four 
special  6w  will  be  obtained  in  the  next  section.  Actually,  a  fifth 
special  case  emerges,  corresponding  to  j  =  -k  and  a  zero  denominator 
k  +  j  ;  6w  for  this  case  can  be  set  to  zero  since  we  still  have  (for 
each  k)  the  unassigned  constant  in  available  for  the  purpose. 


We  start  by  rewriting  (5*0  as 

6w  =  2q "2  ( 6e  sin  v  +  eq_1  6M  cos  v)  +  •£■  eq-2  ( 6e  sin  2v 

+  eq-1  6M  cos  2v)  +  ■§•  e2q”3  <$m  +  q"1  sl  ,  (69) 


where  6e  ,  6M  and  <5L  are  available  from  the  integrals  of  (48),  (50) 
and  (51).  The  integrals  for  6e  and  6M  combine  in  a  very  natural  way 
and  we  eventually  get,  changing  the  interpretation  of  j  (as  in  the 
analysis  for  6r  )  so  that  we  use  the  same  Sj  in  each  term-. 


Sw 


iq  2  A£k 


3e‘ 


'l  -  £  +  k 
,  k  +  j  +2 


+ 


1  +  £  -  k 

k  +  j 


B£ ,  j  +2 


+ 


2e 


1  +  £  +  2k 
.  k  +  j  +  2 


+ 


6 


1  -  £  +  k 
k  +  j  +  1 


+ 


„  1  +  £  -  2k 
5  k  +  j 


+ 


2 


1  +  £  -  k 
k  +  j  -  1 


+ 


L? 1  •  i  i  * 

1  K  *  J  *  2 


+ 


8 


1  +  £  +  2k 
k  +  j  +  1 


„  4(1  -  21)  +  e2( 5  -  £) 
k  +  j 


+ 


8 


1  +  £  -  2k 
k  +  j  -  1 


+ 


e2 


1  +  £  -  k- 

rrrrj 


+ 


1  +  £  +  k  .  1  +  £  +  2k 

f  k  +  j  +  1  6  k  +  j 


♦  6 


1  -  £  -  k 
k  +  j  -  1 


+ 


1  +  l  -  2k  1  n 
k  +  j  -2 


♦ 


3e2 


1  +  £  +  k 
.  k  +  j 


+ 


1  -  £  -  k 
k  +  j  -  2 


B£,  j-2  J  SJ 


(70) 


Though  the  algebra  is  tedious,  we  can  eliminate  j+2  and 
B £  j_2  via  appropriate  versions  of  (29).  If  we  express  ine  result  as 


6w 


•to'?  »tk  l  (Vj.,  B 


1  “t.J.t  *  Vj,0  ®Hj  *  VJ,-1  SJ 


(71) 


the  formulae  for  V,  1  ,  V..  q  and  V,  _1  are  initially  very  complicated. 
They  can  be  greatly  simpl fried,  however;  for  V,  g  this  was  done  by  a 
technique  akin  to  partial  fractions.  The  resulting  expressions  are 

1  6  ^3^2 


■J.1 


'j.O 


-  8 


2e(£  +  j) 
'i  +  2k  +  1 


k  +  J  +  2 


k  +  j  +  1 


k  ♦  j 


k  +  j  -  1 


(72) 


k  +  j  +  1 
-  2e2 


2£  -  1  +  £  -  2k  +  1 


k  +  j 


k  ♦  j  -  1  J 


£  +  k  ♦  1  _  2(  £  ♦  1  ) 

k + j  +  2  k+j 


♦  *-k+1 

k  +  j  -  2 


(73) 
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and 


VJ,-.  ■  2eU  -  J>(irrjTT  *  FT 


j  k ♦ J -  1  k ♦ j  -2 J  • 


As  a  result  of  this  considerable  simplification,  Vj  1  B£,j+1  an' 
V-  _i  •_1  are  now  in  a  form  suitable  for  the  elimination  of  j  +  -| 
ancl  B2,  i-i  in  favour  of  B£,  and  B^.ij  ,  via  an  application  of  two 
recurrence  relations,  (30)  ana  a  similar  one.  Thus,  if  we  now  write 


we  get 


6w  =  TA£k  l  (W£>0  B^j  +  Wjlj_1  bjl--|  , j  )Sj  , 

-  2fi_US_H  .  2  i_i_L  +  1  lJLL±\ 
1,0  ^(k  4  j  +  2  k  +  j  k  +  j  -  2  J 


W£,-1  =  "  2(£  -  1} 


k+j+2  k+j  +  1  k+j 


k+j-1  k  +  j  -  2 


The  final  result  we  require  follows  from  the  substitution  of  (76) 
and  (77)  into  (75).  Writing  Sj  in  full,  we  get 

4"lk  ”  #lk  j  (k  *  j  •  2)  (k  *  j)(k  *  j  -  2)  ([2<1  *  '>  "  k<l<  *  Btj 


.  j  j  -  1)  Bt-1,j}  sln  lku-  *  Jv) 


Equation  (78)  is  the  general  formula  for  6w  due  to  uj£  .  As 
with  (61)  and  (68),  for  6r  and  6b  respectively,  it  applies  for  all 
l  Z  1  ;  like  (68)  but  unlike  (61),  on  the  other  hand,  values  of  |j  |  up 
to  £  -  1  (not  just  £  -  2  )  are  required  to  cover  all  the  non-zero 
terms.  For  each  k  ,  zero  denominators  exist  for  five  different  values 
of  j  :  for  four  of  these  values  (  j  -  -k  ±  1  and  j  -  -k  ±  2  ) ,  spe¬ 
cial  formulae  are  required,  in  place  of  (78),  as  already  noted;  only 
for  the  fifth  value  (  j  -  -k  )  can  a  term  (for  each  k  )  be  actually 
excluded.  It  should  be  noted  that  one  specific  null  term  arises  for 
each  even  value  of  £  .  Thus,  for  k  ■  2  and  j  *  £  -  1  ,  we  see  from 
(78)  that  the  coefficient  of  B^j  is  identically  zero  (independently  of 
£  ),  but  B£_1  1  is  itself  zerowhen  j  -  £  -  1  ,  so  this  specific  term 
of  6w£  2  always  vanishes.  On  this  basis  we  can  derive  the  total 
number  of  terms  required  to  express  6w  for  a  given  value  of  £  , 
counting  in  the  terms  derived  in  the  next  section.  The  result  (derived 
in  the  Report)  is  £2  when  £  is  odd  and  £2  -  1  when  £  is  even. 
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THE  SPECIAL  CASES,  AND  INTEGRATION  CONSTANTS 

The  main  results  in  this  section  are  the  formulae  required  to 
supplement  (78),  the  general  formula  for  6w  .  These  formulae,  covering 
the  cases  j  =  -k  +  1  and  -k  ±  2  ,  are  forced  by  the  'constants’  for 
6e  and  6M  ,  which  are  determined  so  that  the  terms  for  j  =  -k  +  1 
can  be  excluded  from  6r  .  We  also  give  the  formulae  for  the  constants 
associated  with  the  main  elements,  i.e.  for  the  quantities  independent 
of  v  that  are  deliberately  introduced  when  the  v-dependent  components 
of  dfi/dv  ,  for  example,  are  integrated;  by  giving  these  formulae  we 
effectively  define  the  mean  elements  underlying  the  theory.  (The  const¬ 
ants  associated  with  J2  perturbations  are  given  in  Ref.  2). 

Mandatory  Constants  for  6a 

We  go  back  to  (34),  the  original  expression  for  6a  due  to  u£  . 

We  can  expand  the  complete  factor  (p/r)^+1  in  terms  of  the  Bg_+2  j 
(cf  the  expansion  via  the  in  (35)).  On  taking  just  the  term'of 

the  expansion  with  j  =  -k  ,  we  isolate  the  constant  term  that  (for  each 
k  ,  and  a  given  )  is  mandated  by  taking  S'  *  a'  . 

The  result  can  be  written  in  the  form  (for  the  'constant'  compo¬ 
nent  of  6a£k  ) 

6a£k(c)  *  -2a(5"2  A*k  BH+2,k  cos  k“’  •  (79) 

Constants  for  6e  and  6M 


We  have  to  derive  the  formulae  for  6elk(c)  and  6M£k(c)  that 
will  legitimize  our  taking  the  terms  in  6r£k  for  j  -  -k  +  1  and 
-k  -  1  to  be  zero.  These  'constants'  will  complete  the  formulae,  for 
6e  and  6M  ,  given  by  the  integrals  of  (48)  and  (50)  respectively. 

We  start  by  noting  that  (61),  the  general  formula  for  6r£k  ,  was 
obtained  by  combining  the  two  different  denominators  from  (59)  and  (60). 

If  we  do  not  combine  the  denominators,  we  can  rewrite  the  formula  as 

-  ’>p  %  i  („ ,  j  —  -  j  — )  ■  (80> 

The  first  denominator  here  is  associated  with  the  I_  summation  of  (57), 
and  if  this  summation  still  applied  for  j  «  -k  +  1  ,  the  result  would 
be  an  infinite  coefficient  of  Bi_1  _k+1  C_k+1  .  We  actually  want  this 
coefficient  to  be  — 4-(£  -  1)p  A^k  , ’since  it  will  then  neutralize  the 
coefficient  that  arises  without  difficulty  from  the  second  denominator 
in  (80).  The  situation  is  reversed  when  j  -  -k  -  1  and  we  want  the 
coefficient  of  B»_1  _k-1  C_k_i  ,  from  the  second  term  of  (80),  to  be 
■£•( l  -  l)p  A^k  (and  not  infinity;  to  neutralize  the  first  term.  What  we 
do,  therefore,  is  to  obtain  the  coefficients  of  C_k+1  and  C_k_.j  that 
would  apply  in  the  absence  of  the  constants  <5eg,i< ( c )  and  6Md.k(c)  •  we 
can  then  derive  the  appropriate  values  of  these  constants  to  cancel 
these  putative  coefficients. 
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The  details  of  this  analysis  (given  in  the  Report)  are  omitted 
here,  and  we  proceed  directly  to  the  results,  which  can  be  written  as 

6e£k(c)  *  "  t  A£k  leB£,k+2  "  ^  +  k  "  4)b£.,i<  +  1  +  2(£  +  2)eB£k 

-  (£  -  k  -  +  eB£jk_2}  cos  kto'  (81) 

and 

6M£k(c)  -  *  e_1  q  A£k  ^eB£, k+2  ~  (Z  +  k  ’  l4)B£,k*1  ’  2keB£k 

+  (£  -  k  -  k-1  -  eB£>k_2}  sin  kw’  .  (82) 

Constants  for  6i  and  SQ 

We  have  to  derive  formulae  for  6 i^k(c)  and  6 £2^^ / c \  to  legiti¬ 
mize  taking  the  terms  for  j  =  -<  *  1  ana  -<  -  1  in  (68),  the  general 

expression  for  Sb^  *  to  be  zero.  The  analysis  is  somewhat  simpler 

than  that  in  the  preceding  section,  in  spite  of  the  complexity  entailed 
by  the  need  to  work  with  both  k  and  <  . 

As  with  6r^k  ,  we  start  by  observing  that  (68)  was  obtained  by 
combining  two  denominators ,  which  appeared  separately  in  equation  (67). 
When  j  ■  -<  +  1  ,  the  second  denominator  becomes  zero  and  no  longer 
operates;  from  the  first  alone  we  get,  as  the  effective  term  in  (68), 
i£  A^k  _K+1  C*  1  .  When  j  =  -tc  -  1  ,  similarly,  the  first 
denominator  in  (67)  does  not  operate,  and  (68)  effectively  reduces  to 
-J-£  A^k  B^  _  .j  C*K_i  .  These  terms  have  to  be  cancelled  by  the  use  of 
<5i£k(c)  6Q£k(c)  »  with  appropriate  k  .  Again  we  leave  the 

details  to  the  Report,  proceeding  directly  to  the  results,  which  are 

6i£k(c)  *  ±A£k  B£k  cos  ku)' 

and 

6fl£k(c)  *  ^cs"2  A£k  B£k  sin  kw’  . 

Forced  Terms  in  6w 

We  now  have,  for  each  u£  ,  only  one  'constant*  at  our  disposal; 
denoted  by  6u>ik(c\  ,  we  shall  determine  it  (in  the  next  sub-section) 
so  as  to  validate  the  nulling  of  the  term  for  j  -  -k  in  the  formula, 
(78),  for  6w^k  .  For  j  -  -k  ±  1  and  -k  ±  2  ,  on  the  other  hand,  we 
are  forced  to  accept  non-null  terms  that  arise,  via  (5*0,  from  (81)  and 
(82);  we  now  have  to  derive  the  formulae  for  these  terms.  For  each  of 
the  four  special  values  of  j  ,  in  principle  we  embark  on  a  procedure 
that  is  similar  to  that  employed  in  the  derivation  of  $e£i<(c)  and 
6M£k(c)  .  though  more  direct.  In  practice,  however,  insteaa  of 
developing  our  four  special  formulae  more  or  less  ab  initio,  we  start 
four  times  from  the  (final)  general  formula  for  6w  ,  (78),  and  modify 
it  each  time  in  the  appropriate  manner,  replacing  'general'  terms  that 
would  be  infinite  by  special  terms  based  on  <5etk ( c )  and  5M£k(c)* 


(83) 

(8U) 
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As  usual,  we  omit  the  details  (supplied  in  the  Report)  and  simply 
quote  the  results.  Our  four  special-case  formulae  can  be  expressed  as 

5wlk,-k+1  =  t  Alk  _  k  +  2)  B£,,k-1 

♦  (1  -  1)  B^_1>k_1}  sin  (v  +  ku>' )  ,  (85) 

6wlk,-k-1  “  “  "S'  Alk  *  k  +  2)  Bl, k+1 

+  (!l-1)  B^_1>k+1}  sin  (v  -  ku')  ,  (86) 

<h  A£k  {3(8.  +  k  +  5)Bj_ik_2 

-  19(1  -  1)B£_1fk_2}  sin  (2v  +  k«' )  (87) 

4T  A£K  |3(l  -  k  *  5>Bljk,j 

-  19(1  -  1)B£_1(k+2}  sin  ( 2v  -  kw')  .  (88) 


It  remains  to  determine  the  constant,  6u>£k(c)  ,  that  legitimizes 
our  taking  the  term  for  j  =  -k  in  (78)  to  be  zero.  We  already  have 

6M£k(c)  .  given  by  (82),  so  we  only  need  to  determine  6L£k(c)  ,  the 

constant  in  6L£k  ,  for  $w£k(c)  t0  be  known  at  once. 

It  turns  out  (with  the  details  in  the  Report)  that 

6Llk(c)  *  *  kc>  Alk  Blk  sin  k“'  •  (89) 

From  (82)  and  the  definitions  of  L  and  ty  ,  it  follows  that 


^lk(c) 

=  "  t  e  1  A£k 

(e  Bl,k+2  ~  +  k  "  ^)Bl,k+1 

+  2keB£k 

4 

(1  -  k  -  J0B£>k_1  -  e  B£>k_2 

}  sin  kw' 

• 

Finally, 

«lk(c)  is  give? 

■;  by  (844),  so  the  formula  for 

6tolk(c) 

is 

6tolk(c) 

-  -  x  e_1  A£k 

ie  Bl,k+2  "  +  k  ”  A)Bl,k+1 

+  2kes-2 

Blk 

4 

(1  -  k  -  *0B£  k-1  -  e  B£>k_2} 

sin  ku' 

• 

RESULTS  FOR 

The  formulae  of  the  two  preceding  sections  are  valid  for  121, 
the  case  1  «  1  being  trivial.  To  exemplify  the  straightforward  use  of 
these  fomulae,  and  the  complementary  formulae  for  the  rates  of  change  of 
the  mean  elements,  results  for  1-1,  2,  3  and  *4  (together  with  an 
analysis  for  the  exceptional  case,  also  trivial,  1-0)  were  given  in 


6wlk,-k+2  = 

and 

6wtk,-k-2  = 

Constants  for  6w 
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the  Report.  Only  the  results  for  £  =  14  were  new,  however,  so  here  we 
confine  ourselves  to  these.  For  convenience  in  expressing  the  formulae, 
we  define  G  =  ttjt  J14  (R/p) 14  and  f 


=  s2 


We  start  with  the  secular  rates  of  change  for  0  and  <n  .  From 
(44)  and  (45)  we  get,  with  £,  -  4  and  k  =  0, 

t  =  480  Gnc( 4  -  7f)(2  ♦  3e2)  (92) 


and 


i>  = 

-120  Gn(8  -  4 Of  +  35f2)(4  +  3e2) 

9 

(93) 

from  which  it  follows  that 

c 

0 

0 

<\J 

1 

II 

•13 

{4(16  -  62f  +  49f2)  +  9e2( 8  -  28f 

*  21  f 2)  }  . 

(94) 

We  avoid  an  explicit  secular  perturbation  in  R  by  modifying  Kepler's 
third  law.  From  equation  (15)  of  Ref.  7  (or  the  formulae  in  the  Report) 
we  require  (for  the  only  non-zero  harmonic) 

rT2  s3 

=  u { 1  +  288Gq3(8  -  4 Of  +  35f 2) } 

• 

(95) 

Next,  we  cover 
a) .  From  (42) - (47) , 

the  long-period  rates  for  all  the 
we  get,  with  £.  =  4  and  k  =  2  , 

elements 

(except 

£  = 

-480  Gneq2f (6  -  7f)  sin  2W  , 

(96) 

•r 

1  = 

480  Gne2cs(6  -  7f)  sin  2w  , 

(97) 

n  * 

-960  Gne2c(3  -  7f)  cos  2w  , 

(98) 

- 

-240  Gnf (6  -  7f)(2  +  5e2)  cos  2oo 

1 

(99) 

fr  - 

-2400  Gne2qf (6  -  7f)  cos  2w 

(100) 

and 

t  - 

-1680  Gne2qf (6  -  7f)  cos  2u  . 

(101) 

Turning  to  the 
number  of  terms  in  Sr 

perturbations  in  the  coordinates,  we  expect  the 
,  6b  and  6w  to  be  11,  11  and  15,  respectively,  from 

the  formulae  given  earlier.  Starting  with  6r  ,  we  note  that  there  are 
terms  for  k  -  4  ,  k  »  2  and  k  »  0  ,  with  values  of  j  ,  a  priori, 

satisfying  | j  J  i  2  ;  but  for  k  ■  2  we  exclude  j  =  -1  ;  and  for 

k  «  0  we  exclude  j  -  ±1  ,  whilst  the  terms  for  j  =  ±2  are  identical. 
Then  (61)  gives,  corresponding  to  the  three  values  of  k  , 

<5r  -  -2Gpf2  {6e2  cos  2(2u  +  v)  +  35e  cos  (4u  +  v) 

+  28(2  +  e2)  cos  4u  +  1  05e  cos  (3u  +  u)  +  70e2  cos  2(u  +  w)}  ,  (102) 

6r  *=  -8Gpf(6  -  7f){2e2  cos  2(u  +  v)  +  1 5e  cos  (2u  +  v) 

+  20(2  +  e2)  cos  2u  -  30e2  cos  2u>}  (103) 
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6r  =  -24Gp(8  -  4  Of  +  35f2){e2  cos  2v  -  3(2  +  e2)  1  . 


(104) 


For  6b  ,  the  effects  are  for  <  =  3  and  <  =  1  ,  the  a-priori 
values  of  j  being  the  seven  with  | j j  <  3  .  For  k  =  3  we  exclude 
j  =  -2  ,  and  for  k  =  1  we  exclude  j  =  0  and  j  =  -2  .  Then  (68) 
gives,  corresponding  to  the  two  values  of  k  , 

6b  =  -4Gcsf  {4e3  sin  3(u  +  v)  +  35e2  sin  (3u  +  2v) 

+  28e( 4  +  e2)  sin  ( 3u  +  v)  +  70(2  +  3e2)  sin  3u 

+  1 40e ( 4  +  e2)  sin  (2u  +  w)  -  140  e3  sin  3w}  (105) 

and 

6b  =  -4Gcs(4  -  7f){4e3  sin  (u  +  3v)  +  45e2  sin  (u  +  2v) 

+  60e(4  +  e2)  sin  (u  +  v)  -  l80e(4  +  e2)  sin  w 

-  20e3  sin  (2v  -  u) }  .  (106) 


For  6w  ,  the  effects  are  again  for  k  =  4,  2  and  0,  with  the  same 
a-priori  j  values  as  for  6b  .  For  k  =  4  ,  all  seven  j  values  yield 
terms,  of  which  five  come  from  the  general  (78);  for  j  =  -2  we  use 

(87)  and  for  j  =  -3  we  use  (85).  For  k  =  2  ,  the  term  with  j  *  -2 

is  excluded,  whilst  for  j  =  3  ,  (78)  gives  an  example  of  a  ’specific¬ 
ally  null'  term;  there  are  non-null  general  terms  for  j  *  2  and 

j  “  1  ;  and  the  terms  for  j  =  0,  -1  and  -3  come  from  (87),  (85)  and 
(86)  respectively.  Finally,  for  k  =  0  the  term  with  j  -  0  is 
excluded;  the  other  terms  come  in  pairs,  being  ’general’  for  j  *  ±3  , 
from  (87)  and  (88)  for  j  =  ±2  ,  and  from  (85)  and  (86)  for  j  *  +1  . 
Corresponding  to  the  three  values  of  k  ,  we  get 

6w  =  -Gf2{4e3  sin  (4u  +  3v)  +  31e2  sin  2(2u  ♦  v)  +  4e(21  +  5e2)  * 

x  sin  ( 4u  +  v)  +  28(3  +  4e2)  sin  4u  +  26e(7  +  e2)  sin  (3u  +  w) 

+  175e2  sin  2(u  +  w)  +  I40e3  sin  (u  +  3u>)}  ,  (107) 

6w  -  4Gf(6  -  7f)  (2e2  sin  2(u  +  v)  +  4e(5  +  2e2)  sin  (2u  +  v) 

+  5(8  -  7e2)  sin  2u  -  80e(5  +  e2)  sin  (u  +  u) 

-  40e3  sin  (v  -  2w) j  (108) 

and 

6w  =  4Ge(8  -  40f  +  35f2) (2e2  sin  3v  _  3e  sin  2v 

-  6(24  +  5e2)  sin  v)  .  (109) 

It  remains  to  cover  the  induced  short-period  terms  from  the  secu¬ 
lar  and  long-period  rates  of  change^  The  effects  induced  by  the  secular 
variation  are  included  by  adding  (0/n)(v  -  M)  ,  (m/n)(v  -  M)  and 
(R/n)(v  -  M)  to  n  ,  (U  and  R  ,  respectively.  Here  T)  and  u  are 
given  by  (92)  and  (94),  whilst  R  is  the  rate  implicit  in  (95),  so  that 

R  -  144  Gnq3(8  -  4 Of  +  35f2)  ;  (110) 
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The  effects  induced  by  the  long-period  variation,  on  the  other  hand,  are 
allowed  for  via  additional  terms  in  the  expressions  for  6r,  6b  and  6w. 
Using  (52)  -  (5*0,  we  find  that  these  additional  terms  are  given  by 


6r  =  480  Gpe( v  -  M)f(6  -  7f)  sin  (u  +  w)  ,  (111) 

6b  =  240  Ge2(v  -  M)cs  {3(4  -  If)  cos  (v  -  w)  -  7f  cos  (u  +  2eu) }  (112) 
and 

6w  =  240  Ge(v  -  M)f(6  -  7f){e  cos  2u  +  4  cos  (u  +  uj)  -  4e  cos  2w|  . 

....  (113) 

EXTENSION  TO  TESSERAL  HARMONICS 

If  the  Earth's  rotation  rate  is  neglected,  so  that  a'  is  still 
an  absolute  constant  of  the  motion,  the  formulae  that  have  been  derived 
require  surprisingly  little  change.  The  main  change  is,  as  one  might 
expect,  the  replacement  of  all  occurrences  of  A£k  and  A£lC  by  quanti¬ 
ties  A£mk  and  A£m<  ,  to  reflect  the  generalization  of  J£  to  J£m  , 
with  the  arguments  of  Cj  and  Sj  ,  in  equations  (1),  replaced  by 
jv  +  ku'  +  m(fi'  -  X£m),  where  O'  «  0  -  v  -  ,  the  sidereal  angle  (v) 

being  supposed  fixed;  further,  we  now  have  to  allow  for  negative  values 
of  k  ,  which  takes  alternate  values  from  -l  to  l  .  (For  amplifica¬ 
tion  of  these  remarks,  with  definitions  of  A£Jnk  and  Ag^  ,  and  a 
discussion  of  two,  possible  definitions  of  the  A£m(i)  functions  that 
generalize  the  A£(i)  ,  see  Appendix  A  of  the  Report). 

By  making  use  of  the  quantities  Ag^  and  Ag,,,^  ,  we  find  no 
difficulty  in  extending  the  theory,  largely  because  the  treatment  of 
(p/r )^+1  ,  via  the  Bgj  ,  goes  through  unchanged.  Thus,  equations  (35), 
(48),  (39)  and  (50),  for  6a  ,  de/dv  ,  diji/dv  and  dM/dv  ,  respectively, 
are  unchanged  apart  from  the  appearance  of  Agm^  in  place  of  A£m  . 
Equation  (38),  for  dn/dv  ,  requires  a  corresponding  change,  such  that 
the  derivative  Aj^  replaces  Aj^  .  This  just  leaves  (49),  for  di/dv  , 
for  which  a  slightly  more  complicated  expression  is  now  required,  to 
reflect  the  fact  that  pc2  is  no  longer  an  invariant.  We  have,  in  fact, 


-  2mpc  Atok  l  By  Sj  , 


(114) 


from  which  we  derive 


■57  ’  s_1(kc  -  m)  Ag^  [  Bgj  Sj  ;  (115) 

in  comparison  with  equation  (49),  we  see  that  the  only  additional  change 
is  the  replacement  of  kc  by  kc  -  m  . 
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Six  of  the  seven  formulae  that  define  6r,  6b  and  6w  completely, 
for  the  zonal  harmonics,  apply  immediately  to  the  zonal  harmonics,  so 
long  as  replaces  A£_m  and  the  trigonometric  argument  includes  the 

term  m(fl'  -  X£,m)  .  These  six  are  (61)  and  (78),  for  the  general  6r  and 
6w  ,  and  (85)  -  (88),  the  four  special  formulae  for  6w  .  In  the  seventh 
formula,  (68)  for  6b  ,  i  Aj,*  must  be  replaced  by  (l  +  m)A£mK  ,  in 
addition  to  the  inclusion  of  the  new  term  in  the  trigonometric  argument. 

The  numbers  of  terms  in  6r,  6b  and  6w,  for  a  given  J$_m  ,  are 
greater  for  m  >  0  than  for  m  =  0  ,  to  reflect  the  distinction  between 
positive  and  negative  k  .  These  numbers  are  otherwise  independent  of 
m  ,  however,  in  consequence  of  which  we  find  that  there  are  2 £2  -  3£,  +  1 

terms  for  6r  ,  2£,2  -  3i,  +  2  terms  for  6b  ,  and  2£2  (for  odd  l  )  or 

2S.2  -  1  (for  even  SL  )  terms  for  6w  . 

CONCLUSION 

The  original  idea  of  Kozai®  -  that  the  short-period  perturbations 
due  to  J2  can  be  more  compactly  expressed  by  amalgamating  6a,  6e,  6w 

and  6M  into  6r  and  5u  -  has  been  carried  to  its  logical  conclusion  by 

eliminating  the  short-period  perturbations  in  all  the  elements  in  favour 
of  effects  in  spherical  coordinates.  The  general  formulae  for  these 
effects  are  given  by  the  summations  in  (61),  (68)  and  (78),  for  the  per¬ 
turbations  in  r,  b  and  w,  respectively.  Terms  that  would  have  a  zero 
denominator  are  excluded  from  these  summations,  as  a  consequence  of  the 
optimal  definition  of  mean  elements,  except  that  replacement  terms  are 
needed  for  the  perturbations  in  w  ;  the  formulae  for  these  are  (85)  - 
(88). 

The  formulae  for  coordinate  perturbations  are  complemented  by 
(i)2)  -  (i)7),  which  are  the  formulae  for  the  rates  of  change  of  the  mean 
elements.  Relatively  speaking,  these  formulae,  which  lead  to  the 
secular  and  long-period  perturbations,  contain  very  few  terms,  but  over 
periods  of  time  longer  than  an  orbital  revolution  the  effects  are  much 
greater  than  those  from  the  short-period  perturbations.  Formulae  to  a 
further  order  of  approximation,  i.e.  to  and  J2  J£  (i  >  2),  have 

been  given  by  Berger  and  Walch9,  and  Kinoshita10,  but  expressions 
appropriate  to  the  mean  elements  used  in  the  present  paper  have  not  yet 
been  derived. 
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